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I. INTRODUCTION 



The excitation spectrum of a two-dimensional px + ipy fermionic superfluid, such as a thin film 
^ ' of ^He-A, includes a gapless Majorana-Weyl fermion which is confined to the boundary by Andreev 

O I reflection. There is also a persistent ground-state boundary current which provides a droplet con- 

taining N particles with angular momentum hN /2. Both of these boundary effects are associated 
5— ( ' with bulk Chern-Simons effective actions. We show that the gapless edge-mode is required for the 

Oh' gauge invariance of the total effective action, but the same is not true of the boundary current. 

^3 ■ 

^ . PACS numbers: 

> ■ 

fi 

Ch ' Read and Green Q have pointed out close parallels between two-dimensional Px + ipy chiral fermionic superfluids 
O (a thin film of ^He-A for example) and the v — \/2 Pfaffian quantum Hall state. The many-body ground-state 
wavefunction of both systems contains a Pfaffian factor, and both support eapless Majorana-Weyl edge modes 0]- In 
addition, both systems can host vortex defects with non-Abelian statistics I^Q. The Px + iPy superfluid also possesses 
an equilibrium edge current that is consistent with an ?i-per-Cooper-pair intrinsic angular momentum 0, ■ This 
current is similar to that induced by a confining potential at the boundary of a Hall droplet, and occurs because chiral 
fermionic superfluids manifest an analogue of the quantum Hall effect, even in the absence of an external magnetic 
_ field % 

The principal differences between the Pfaffian Hall state and the chiral superffuid are that, in addition to the 
Pfaffian, the Hall-state wavefunction includes a symmetric v ~ \/2 Laughlin factor; also the Hall system is spin 
CO polarized, while the superfluid retains two active spin components. A less signiflcant difference is that the Hall fluid is 
incompressible while a neutral superfluid like '^He-A has a gapless acoustic mode. This acoustic mode will, however, 
be gapped in a charged superfluid state, such as that believed to exist in the layered Sr2Ru04 superconductors PI.ITol|. 

In the quantum Hall effect, the lack of low energy bulk degrees of freedom leads to the dynamics of the boundary 
and the bulk being intimatel y c onnected. Currents flowing in the the bulk serve to soak up the anomalies of the 
edge- mode conservation laws and, more profoundly, the bulk many-body wavefunction can be constructed from 



l y C O 

,0, ^ ^, ^ ^ 

the conformal blocks of the edge conformal field theory (3j. The existence of the conformally invariant Majorana edge 



I 

Q , mode suggests that something similar is true for the superfluid. 

O ' Because of the interest of quantum systems with non-Abelian statistics, both for their intrinsic appeal as exotic 
physics and for their potential application to quantum computing .12^ .13i | , there is reason to search for the simplest 
possible description of their dynamics. Thus we seek an effective action which captures all the essential low-energy 
' degrees of freedom of the system. The construction of such an effective action for the fermionic Pfaffian Hall state is, 
however, rather complicated |l4j . It is therefore worth examining other systems with similar properties. 

One such system is the bosonic Pfaffian quantum Hall state, which can perhaps be realized in rotating Bose 
condensates. The bosonic Pfaffian state has a ground-state wavefunction consisting of the product of the antisymmetric 
Pfaffian and an antisymmetric = 1 Laughlin factor, making it even under particle interchange. In this case i4,il the 
effective action is an SU{2) Chern-Simons term at level k = 2. This means that the edge states form representations 
of an SU{2)2 current algebra. Indeed, the space of low energy states is spanned by wavefunctions consisting of 
polynomials in the complex coordinates Zi that vanish when any three of the Zi coincide. The generating function for 
the number of such polynomials is a character of the SU{2)2 algebra [T^ . 

The px + ipy superfluid is also relatively easy to analyze because its effective action has been computed by standard 
gradient expansion methods ^3 ■ It the purpose of this paper to explore some of the mathematical and physical 
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aspects of the resulting expression. We will see that the response of the fluid to a non-Abelian gauge field which 
couples to the spin degree of freedom is governed by a Chern-Simons term of conventional form, but with a coefficient 
corresponding to level k — ^. Since invariance under "large" gauge transformations demands that the level k be an 
integer, this means that the spin action cannot be gauge invariant on its own. Gauge invariance is restored only when 
we take into account the role of the spin-triplet order parameter in spontaneously breaking of the SU{2) symmetry 
down to C/(l), and then by the subsequent absorption of the residual U{1) gauge dependence by the C/(l) axial 
anomaly of the chiral edge-states. The topological origin of the edge modes is therefore illuminated. The response 
to an Abelian particle-number gauge field is governed by an action which looks superficially Chern-Simons like, but 
is deficient in that it lacks the terms with time derivatives. In this case, gauge invariance is immediately manifest 
once we include the U{1) Goldstone mode in the effective action In particular, no boundary degree of freedom 
is required. Although this means that there is no explicit coupling of the Abelian gauge field to the edge modes, the 
Abelian part of the effective action does describe the Hall- like edge current, which therefore has a different origin 
from the edge modes. It also reveals an essential difference between the real Hall effect and its field-free analogue. In 
the Hall effect the current is proportional to the applied electric field. In the superfluid, the absence of the A in the 
action means that the induced current is proportional not to the electric field, but to the gradient of the fiuid density. 
The "Hall current" should therefore be thought of as a two-dimensional analogue of the Mermin-Muzikar current |y] , 
which is due to the intrinsic angular momentum of the fiuid. 

There have been several recent papers discussing the effective actions for chiral superfluids and superconductors 
[isL , but our point of view differs from these in its interpretation of the "Hall effect" , the magnitude and origin of 
the edge-current, and in its emphasis on the role of gauge invariance in establishing the bulk/edge connexion. 

In the next section we will discuss the Bogoliubov action for the planar + ipy spin-triplet superfluid. In section 
mil we will find the eigenfunctions of the corresponding Bogoliubov-de Gennes Hamiltonian in rigid walled containers, 
and use them to compute the edge-mode spectrum and the magnitude of the persistent edge-currents. In section four 
we will describe the effective actions governing the response of the fluid to external gauge fields that couple to particle 
number and to spin. We then show how the existence of the persistent boundary current and gapless edge-modes can 
be deduced from the bulk effective action. 

II. BOGOLIUBOV-DE GENNES OPERATOR 

The fermionic part of the action describing a BCS superconductor can be written in Nambu two-component for- 
malism as 

S= f d^xdt (]--^\idt - i?)* ) . (1) 
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Here, 

* = 



*^ = [V'L V-a], (2) 



are Grassmann fields with a spinor index a =t; i, and H is the Bogoliubov-de Gennes Hamiltonian 



H 



h A 



(3) 



The entries in H are single-particle operators acting on the tensor product of the position and spinor spaces. The 
hermiticity of H requires that the single-particle Hamiltonian h be Hermitian, and Fermi statistics requires that the 
gap function A be skew-symmetric, in this combined space. 

For a Galilean invariant fluid of particles of mass m, the Hamiltonian is 

where / is the identity operator in spin space and e / the Fermi energy. For a superconductor h can be a more general 
function, e(p)/, of p = — iV. 

The gap function will usually be a dynamical field, and the total action will contain additional terms that serve to 
determine its value via a gap equation. We are, however, interested only in the response of the fermions to prescribed 
changes in the gap function A, so we will not need to make these terms explicit. Further, we are primarily interested 
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in topological effects, so the precise form of A is not significant as long as it produces the required symmetry breaking. 
For & + ipy, spin-triplet superfluid, such as ■^He-A, we may take 



,i4./2 



(5) 



Here, { , } denotes an anticommutator, A is the magnitude of the induced gap in the quasiparticle spectrum, and 
$ is the overall phase of the order parameter. The spin part S is a symmetric 2x2 matrix with entries 



= • d)cT2)a/3, 



(6) 



where d is a unit vector. The orbital part of the gap function is contained in the operator P, which we will take to be 



P = -i{px + iPv) = -{dx + idy), 



(7) 



corresponding to Cooper pairs with their 1=1 angular momentum vector 1 directed in the +z direction, i.e. perpen- 
dicular to the plane of the fluid. This orientation for 1 ensures that the entire Fcrmi-siirfacc is gapped. We will fix 
1 = z throughout this paper. We will also take the magnitude of the gap A to be a constant. Although this parameter 
should be determined self-consistently through a gap equation, its magnitude serves only to provide an upper limit 
for what we mean by "low-energy" degrees of freedom, and so any variation has no role in the following discussion. 

When E and $ are constants, the operator ordering of the spin, phase, and orbital parts of A is unimportant. When 
these quantities vary in space, however, we need the anticommutator of S and P, and the symmetric distribution of 
the overall phase e** about it, to ensure the antisymmetry of A. 

In the sequel, our calculations will be strictly 2-1-1 dimensional. They therefore apply to a single stratum of a 
layered superconductor, or, for a thin fllm of "'He-A, to the case when only a single transverse momentum mode lies 
below the Fermi surface. When n transverse modes are occupied but the order-parameter remains independent of z, 
our eflfective actions should be multiplied by n. 

III. EDGE MODES AND BOUNDARY CURRENTS 

We now suppose the superfluid to be confined by rigid walls at which the wavefunction is required to vanish. We 
also assume, for the duration of this section, that the spin vector d lies in the y direction, so that Y, = il and the two 
spin components decouple. 



A. Rectangular Geometry 



If we substitute 



(8) 



with constant a, b, into the Bogoliubov equation, = E^, the eigenvalue condition reduces to 



E 



e{k) 

(^)e-^A -e{k) 

Here 9 is the polar angle such that kx = k cos 9, ky = k sin 9. The plane- wave eigenstates are therefore 



(9) 



i(T3e/2_ 



2^E{E + A) 



E + e{k)+A' 
E - e{k) + A 



(10) 



with E = +^e'^{k) + A'^{k'^/k'j), and 



2^\E\{\E\-A) 



\E\-A-e{k) 
\E\-A + e{k) 



ikxX-\-ikyy 



(11) 
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r(k,E) 



FIG. 1: Geometry of specular and Andreev reflection from the boundary. 



with E = —^e'^{k) + A'^{k'^/k^). We will usually be interested in states close to the Fermi surface where we can 
approximate the energy as 



(12) 



where w/ is the Fermi- velocity. We have considered only one spin component. There are really two sets of such 
solutions, one for spin up and one for spin down. 

If we expand the quantized field operator \l/ in terms of the plane wave states 



(13) 



then, in order for the upper and lower components to reconstruct t/j and V'^ respectively, the operators a^j^k must obey 
the reality condition aE,u — oL^ -k- 



Edge states 



We now investigate the effect of a boundary. Suppose, as shown in figure ^ that there is a wall at y 
fluid lying above it. We seek solutions to the Bogoliubov equation in the form 



b+iy) 



^ikfX cos 9 -\-ikfy sin ^ 



a~iy) 
b-iy) 



^ikfX cos 9—ikfy sin 6 



0, with the 



(14) 



where now a± and b± are allowed to vary with y — but slowly on the scale of fc/. We impose the boundary condition 
that ^ = at y = 0. The resultant equations for a± and b± coincide with those derived in the appendix for the 
solutions of a one-dimensional Dirac hamiltonian. The vanishing of the wavefunction at the wall becomes (because 
of the minus sign before C_) the continuity of the Dirac wavefunction at x = 0. The incoming particle therefore 
sees the reflecting boundary only as an abrupt change in the phase of the orbital part of the order parameter, which 
jumps from —6 to +9. The transmitted Dirac wave with amplitude t{k,E) corresponds to those particles that have 
been specularly reflected at boundary, and the backscattered Dirac wave with amplitude r{k, E) corresponds to those 
particles that have been Andreev reflected, and so retrace their path. The only significant differences between the 
two dimensional geometry and the one-dimensional problem solved in the appendix arc that the co-ordinate "a;" in 
the appendix is the distance along the classical trajectory, and that there we have set the Fermi velocity Vf to unity. 
Thus, expressions appearing in the appendix, such as e"'''^' must here be replaced by e"'^'^'/^''^'''"^^)^. 

On making the appropriate translations from the results in the appendix, we find that we have a bound state 



Xo(fcx) 




ikrx cos 6 



3in(fc/y sir 



with energy 



So(fcx) 



-A cos 61 



-A{k^/kf). 



(15) 



(16) 



The associated edge excitations are therefore chiral, or Weyl, fermions, having group velocity dE/dk^ only in the — x 
direction. This motion is in the opposite sense to the anticlockwise orbital angular momentum of the Cooper pairs. 
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These edge modes have a topological character. All that is required for them to exist is that the fluid density fall 
to zero at the boundary 0. Indeed the edge- mode spectrum 1)16(1 for our rigid wall boundary coincides with that 
obtained by Read and Green Q in the opposite extreme of a very soft wall. 

The contribution of the edge-modes to the field operator is 

X]ao,fc,Xo,/c,(a;), (17) 

and, in order for the upper and lower components to reconstruct -0 and tp^ respectively, we must have ao,fc^ = ■ 
The edge-mode field 

is therefore real, = ^q{x), or Majorana. 

It is difficult for anything interact with a Majorana- Weyl particle. A Weyl particle can interact only via currents, 
and no current operator can be constructed from a single Majorana field. In our superfiuid there are, however, two 
Majorana- Weyl fields: one for spin up and one for spin down. Out of these two real fields we can construct one 
complex Weyl field 

*c = -^(^'To + i'J'io), 

^1 - ^(*To-^*io), (19) 
and from this we can construct a unique U{1) current operator 

^I'^c = ^ (*TO^'io - ^'io^'To) = ^*of^2^'o- (20) 
The edge modes may therefore interact with the (72 component of a spin-coupled gauge field. 



2. Boundary current 

The doubling of the number of degrees of freedom in the Bogoliubov formalism requires us, when computing ground- 
state expectation values of an operator, to sum the contributions of all occupied (negative energy) states, but then 
divide by two to prevent over-counting — this division by two being equivalent to imposing the reality condition on 
the aE,]i- The mass current carried by the state 

^ ^ ^ik fX cos 9-\-ik fy sin 9 (21) 

where a, b are slowly varying compared to kf is therefore 

j = i(|a|2 + |6|2)(/c^cos6l,fc/ sin 61). (22) 

The edge modes with fc^ > have negative energy, and so are occupied in the ground state where they make a 
positive contribution to the boundary momentum density. For a Galilean invariant system, the momentum density is 
also the mass current. The edge modes therefore tend to produce a boundary current that fiows in the same sense as 
the Cooper-pair rotation. The bound states, however, are not the only contribution to this boundary current. The 
balance is provided by the phase-shifted scattering states, which (as they do the theory of fractional charge [23ll2^) 
partially cancel the contribution of a bound state when it has negative energy and is occupied, and partially make up 
for the absence of its contribution when it has positive energy and is unoccupied. 

The current in the direction perpendicular to the wall will cancel between the incoming and outgoing waves j2^, 
but the X components will add. Using the results from the appendix, we find that the net mass-current running near 
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the boundary is (for a single spin component) 



^^^^ ^ / . X ^^^^^^ f / . _ r ^^^^^ 



"'0 



27r \2tt J ' I, 2tt \2tt 2 



d(kf cos 6) e 

= 2 / L kf cose 

Jo 2^ 27r ^ 

fc? 1 

where p — k'j/An is the number density per spin component. If this current flows at the edge of a disc-shaped region 
of radius R, it provides angular momentum 

L = 2nR' (9 z = f z. (24) 

This result agrees with that of Kita 7] , and Volovik [2nl | and is exactly the angular momentum we would expect of a 
fluid of tightly bound Cooper pairs, each pair having orbital angular momentum of hz. That the same result holds in 
the weak coupling limit, where only the particles near the Fermi-surface are affected by the pairing, is more surprising, 
and it is often referred to as the "angular momentum paradox" (see Kita for a recent review of this). 

Our result for the boundary current differs from that of Furusaki et al. jioj , because they consider only the the 
bound state contributions. The boundary current arising from the bound states alone is (again for a single spin 
component) 

1 rUjk.cose) kj^ejrn^ (25) 

2 io 27r 87r 47r ^ ' 

After multiplying by two to take into account the two spin components, this coincides with equation (2.8) in Furusaki 
et al. p^ . and differs from the actual persistent boundary current by a factor of two. 

We have computed the boundary current only in the weak coupling limit, but the result that it precisely accounts 
for the fi.-per-pair angular momentum should remain valid even as the coupling is increased. This is because, provided 
the ground state evolves adiabatically, its total angular momentum cannot change as we alter parameters in the 
Hamiltonian. Adiabatic evolution may fail if there is spectral flow, and such flow does occur and change the angular 
momentum when we compute the angular momentum of some vortex configurations in bulk ^He-A p]| . and also when 
we consider the dynamical intrinsic angular momentum density |22j| . but in the present case spectral flow can only 
occur through the gapless edge states and the Eo{kx) ~ —Eo{-~kx) Majorana symmetry prevents the spectrum moving 
en-masse with respect to the chemical potential. For the same reason, the total boundary current should also not be 
affected by local deviations of the order parameter away from its bulk form. The "fractional charge" interpretation 
123 of the current ensures this, because the total fractional charge depends only on the asymptotics of the 
order parameter, and is unaffected by local variations. 

B. Disc geometry 

When the fluid has disc geometry it is convenient to use polar coordinates r, 0, and to work with angular momentum 
eigenstates oc e*'*. For example, we have 

- V^Ji{kr)e'^^ = k^Ji{kr)e'^'', (26) 
where Ji{kr) is a Bessel function. In polar coordinates, the orbital part of the gap operator P = — {dx +idy) becomes 

and has a particularly simple action on the Laplace eigenfunctions: 

Pe'^^Mkr) = fce*('+i)«Ji+i(fcr). (28) 

The adjoint operator 
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similarly acts to reduce the angular momentum eigenvalue 
We look for eigenstates in the form 



bJi{kr)e^^'> 



The Bogoliubov equation, H'i> — E'i, then reduces to 



e(fc) {k/kf)A 
(fc/fc/)A ~e{k) 





a 


= E 


a 










b 


b 



The eigenstates are therefore 
where the energy eigenvalue is 



+ e + A)e'('+i)«J,+i(fcr) 
{E - e + A)e'^^ Ji{kr) 



E 



and 



with energy 



2^\E\{\E\-A) 



(1^1 - A + e|£:|)e*'^J,(fcr) 



E = -^e-^{k) + {k/kfY/S?. 
As before, we are interested in momenta near the Fermi-surface where the energy can be approximated by 



E = ±^lv){k-kfY + A? 



(30) 



(31) 



(32) 



(33) 



(34) 



(35) 



(36) 



(37) 



When we confine the fluid by imposing rigid wall boundary conditions at r = i?, the eigenstates will be a linear 
combination 



^i = C+ 



a+Jz+i(fc+r)e'('+i)« 
fe+Ji(fc+r)e''« 



a_J,+i(fc_r)e^('+i)^ 
b^Ji{k^r)(^^^ 



of unconfined states with slightly different momentum k± = kf ± k, but a common energy 



The coefficients a± and b± are given by 



E 


= ± 


a± 









E±vfk + A 
E^vfk + A 



(38) 



(39) 



(40) 



To examine the consequences of the condition that ^' = at r = i? we use the WKB approximation for the Bessel 
function 



Ji{kr) 



TTkx{r) 



sm{kx{r) - W{r) - tt/A), r':$>b. 



(41) 



Here x{r) and 9{r) are functions of r, defined in terms of the semi-classical impact parameter, b — l/k, by a; = j'sin^ 
and b — rcos9. As illustrated in figure [3 the parameter x has the physical interpretation of being the distance along 
the straight-line semi-classical trajectory. This approximation is quite accurate once r exceeds b by more than a few 
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FIG. 2: The geometry of the WKB approximation to the Bessel function. 

percent, and is therefore reliable except for a few large values of I which correspond to classical trajectories grazing 
the boundary. Using the WKB approximation and the explicit form of the coefficients a± and b± we we end up with 
exactly the same equations for the bound state and S matrix that we found in the planar boundary case. 

There is some advantage to working with a circular container, however. With a finite length boundary, the set of 
edge modes becomes discrete — being labeled by the integer I. The angle 9 between the semi-classical trajectory and 
the boundary is 9{R), and I = kfRcos9{R). In terms of I, therefore, the bound state has energy 

Eo{l) « - (r^) A, (42) 

where Imax = fc/i?. The WKB approximation is not quite good quite good enough to distinguish between I and I + ^ 
in this expression, but on general grounds we know that if (u, v)^ is an eigenstate of the Bogoliubov Hamiltonian 
with energy i?, then (w*, u*)^ is also an eigenstate with energy —E. Under this transformation '^i and so 

the correct equation must be 

^.(0 . - (ifiZH) A. (43, 

There is therefore no exact zero-mode. If, however, the bulk fluid were to contain an odd number of vortices, and thus 
the phase of the order parameter wind an odd number of times as we encircle the boundary, then the Vs appearing 
in the upper and lower components of \E' would differ by an even number. In this case there will be a zero energy 
edge state. Since each vortex has a zero- mode in its core 0, there will be an unpaired zero-energy core mode which 
can pair with the zero-energy edge state and so preserve the even dimension of the total Bogoliubov-particle Hilbert 
space. 

This dependence of the edge-mode spectrum on the number of the bulk vortex excitations is reminiscent of what 
happens in the quantum Hall effect. There the Hilbert space sector of the edge conformal field theory also depends 
on the number and type of vortex quasi-particles that are present in the bulk. 

IV. EFFECTIVE ACTIONS 

We now discuss the effective action governing the low-energy dynamics of the superfluid. We will obtain the action 
by examining the response of the fluid to external gauge fields that couple to the particle number and spin of the 
fluid. 



A. Particle-number symmetry 

We begin by gauging the C/(l) symmetry corresponding to particle-number conservation. We will again hold the 
d vector fixed in the y direction so we can treat each spin component separately. We minimally couple the particle- 
number current to an Abelian gauge field (Aq, A), where Aq is the time component and A = (Ai, A'l) are the in-plane 
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components of the externally imposed field. This requires the replacement 

idt-H ^ idt - H{A, $) 

where 



F(A<f) = 



^0, 



-JA e-*/2pte-*/2, i(v + zA)2 + ylo 



The resulting action 



is then invariant under the local gauge transformation 



provided we simultaneously transform 



(44) 



(45) 



(46) 



(47) 



A 

^0 



$ + 20, 
A + V(/., 
Ao + dtcf). 



(48) 



Because the Grassmann measure in the path integral is left invariant by the transformation H47|l , the effective action 

iS^^'^iA,^) = In exp[iS'(A, 5'^)] 



■InDet {idt-H{A,^)) 



(49) 



will be invariant under the transformation (|48|l . If we compute H49|l to second order in the gauge field and its gradients, 
we find [^Il8||. for a single spin component, 



5e"r(A*) 



d^xdt 

-<7xy I 



P0_ 

2m 
fd<S>/2 



1 /9$/2 



Ao 



- Ao (V X A)^ - po 



(V$/2- A) 

/a$/2 



V dt 



-Ar 



(50) 



Here, (T^y — g^, the parameter Cg is the speed of sound, and po the equilibrium number density. The coefficient a^y 
is topologically stable, but the other quantities depend on the details of the fluid 8]. For a 2+1 dimensional Galilean 
invariant system of particles with mass m we have po = and = i^//2. The action H5U|) is manifestly invariant 
under the transformation (|48|) . 

The term with coefficient <7xy contains a Ghern-Simons-like part 



(fix dteoijAo diAj . 



(51) 



This is not a complete Chern-Simons action, however, because there is no SiojAi dtAj term. It does, nonetheless, 
imply the existence of a Hall-like response to the external field. We find for the particle-number current 



6A m 

PqVs + crxy(z X V) 



(V'I>/2 - eA)) + V) 



a$/2 

dt 



Ao 



\ dt 



(52) 
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Although the term with axy contains a gradient of Aq, it is not equal to <Txy{E x z), where E = V^o — A, as it would 
be in the Hall effect. (Observe that V<i>/2 cannot be A in disguise, because the former is necessarily curl- free.) We 
note, however, that the combination (9t$/2 — Aq) occurs in the the expression for the density 



r cnum 
' Mo 

Consequently, it seems preferable to write 



mci 



(JxyV X A. 



(53) 



PoVs - — (z X V)(p - CFxyBz), 



(54) 



where — x A)z, and so recognize that the "Hall" current depends on the external field primarily through its 
effect in modifying the density. The natural analogy is then with the bound current j bound = V x M in a magnet with 
varying magnetization M. In the superfluid, the magnetic-moment density M is replaced by the intrinsic angular 
momentum density ■^fi{p — axyBz)i. The cTxyBz correction is presumably present because a diamagnetic response to 
the external field will reduce the kinetic angular momentum of a Cooper pair even while leaving its canonical angular 
momentum fixed at h. In the absence of an external field, we therefore have a mass flux 



(55) 



and this is a planar analogue of the Mcrmin-Muzikar current |^] . When p goes to zero slowly at a boundary we can use 
this Mermin-Muzikar expression to compute the equilibrium boundary-current. The resulting boundary momentum 
density coincides with that we computed for a rigid wall in section llTTl and again provides an ft-per-Cooper pair angular 
momentum. This agreement between the rigid-boundary calculation and the gradient expansion when expressed in 
terms of the density, coupled with the physical explanation in terms of the intrinsic angular momentum of the fluid, 
leads us to conjecture that the total mass-current is always proportional to the change in density, rather than to the 
external force that causes the change. 



B. Spin- rotation symmetry 



The fields in 



S = J Sxdti^^idt-H)-^ 



(56) 



can also be coupled to an SU{2) gauge field which acts on the spin indices. To do this we replace the derivatives in 
S by covariant derivatives 



"t" An, , 



where 



An icJaAn 



^p, U- -fj, 

is an externally imposed gauge field. Under a local gauge transformation the Fermi fields transform as 



where U € SU{2). The covariant derivatives transform as 

d,, + A,, ^ U-\d^ + A^)U = 9^ + {U-^A^U + U-^d^U). 



(57) 
(58) 

(59) 

(60) 



When considering how the transformation act in the — h?" entry in , we need to recognize that derivative operators 
appearing there are the transpose 9j = —9^ of those in h, and so the covariant derivatives will be also be the transpose 
+Al = -df, + Al. Using U* = {U-^Y = (C^^)"\ we have 



= -d^ + {u-^A,,u + u-^d^uf 
= [d,, + {u-U^u + u-^d^u)f, 



(61) 
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and the the effect of the transformation is consistent for both entries: 

^ = C/-U^C/ + U-^d^U. (62) 

Note that {A^)^ — A^^ . For the off diagonal terms we have 

U''^Y.U* = U-^ {i{a ■ d)(72) U* = {iU-\a ■ d)Ua2) ■ (63) 

The net result is that the gauged action is invariant under the transformation (|59l) . provided we simultaneously 
transform 

(d-cr) i7(cr-d)[/-l 

A^ ^ A]^'" - UA^U-^ + Udf,U-\ (64) 

Following Volovik and Yakovenko [l3|, we now compute the effective action with the d vector fixed to lie in the y 
direction, and find a Chern-Simons term 



Sir id -y^A)^^Jtr (^AdA+lA''^ . 



(65) 



Here we use differential-form notation in which A = iaaAl^dx'^ is a matrix valued 1-form. For any compact simple 
gauge-group G, the Chern-Simons action is defined to be 

C{A) = ^ y" (^AdA + ^A^^ , (66) 

where A = iXaA^^dx^ is a Lie (G')-algebra- valued 1-form, and, as is customary, we normalize the trace and the 
Hermitian generators Xa by tr (AaAb) = 26ab- When such a term is to appear in a functional integral. 



d[A]e'''^'^^>, (67) 

then coefficient k must be quantized. This is because under a gauge transformation 

A^A3 = g-'Ag + g-'dg (68) 

we have 

C{A) -> CiA^) = C{A) - 7^ / tr [ig-'dgf] - ^ / tr (dgg-U). (69) 

Here, is the region occupied by the fluid and Oft its boundary. The last term is zero when fl is closed manifold, or 
if we restrict the gauge transformation to those constant on the boundary. The second term, 

W(5) = / tr [{g"'dgf] , (70) 

has no reason to vanish, however. It is the pull-back to fl of an element of, H^^{G, Z), the third dc-Rham cohomology 
group of G. It can be non-zero whenever the gauge transformation g{x) maps £7 into a homologically non-trivial 3- 
manifold in G. In particular, when VL is and the group G is SU{2), then 

-|- / tr [{g-'dg)'] = 27rn, (71) 

where n is the degree of the map from SU{2) ~ . In this case k has to be an integer so that the gauge 

ambiguity in G{A) is 27rfcn and exp(zA:C(^)) is well-defined. 

The Chern-Simons action found by Volovik and Yakovenko has a coefficient corresponding to k — 1/2, and so 
violates the quantization condition on k. It cannot be gauge invariant on its own. The complete effective action will 
be gauge invariant, of course, but we must include additional degrees of freedom to make this manifest. One of these 
is the direction of the vector d, which we must therefore allow to vary. We parameterize d in terms of a group element 
V e SU{2) by setting 

{d- a) = Va2V-^. (72) 
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Then, 

Sirid.A) = 5:r(y,^^) ^ IciA""). (73) 

This expression is clearly invariant under the simultaneous replacement A — > and V — > U~^V . This is good, 
but not perfect. The problem is that V is not unique. The vector d is more correctly paramaterized by elements of 
the coset S'C/(2)/C/(l), since we can replace V by without changing d. This substitution does affect C{A^), 

however. Compensating for the effects of e*°'^"^ requires yet another degree of freedom. We will write this as = e^"'^^. 
A completely gauge invariant action is then 

5:r(d,Ax) = ic(^^^). (74) 

This is manifestly invariant under the simultaneous replacement 

V -> U-^Ve"""^, 

X^X-(I>- (75) 

What is the physical interpretation of this extra field x? This question is easiest to answer if we set V — I o\\ the 
boundary, so that d = y there. Then, using the Polyakov-Wiegmann identity (|69|l we can write 

STi^^Ax) = ic(^^^) = ic(A)-^ / tr [{V-'dVf] 



-— / tr [dWW-^A] 
Stt Jan 



(76) 



The second term 



^m^) = ^^tr [[V-^dVf], (77) 

is precisely the Hopf index found by Volovik and Yakovenko. On a closed manifold, or when d is fixed at the boundary, 
it is equal to titt where n labels the homotopy class of the map d : 5^ — s- S'^. The Berry phase provided by this term 
makes a skyrmion soliton in the d field into a fcrmion. The third term, 

- — / tr {dWW~^A] , (78) 

is equal to 

~i^ldn ^'^^'""-^^ = i / '^'"^^ ^-^"^"^ " ' ^^^^ 

and represents the coupling of the x field current to the cr2 component of A. This interaction takes place only on the 
boundary, which we have taken to be the x axis as in section UlI Al This suggests that the x field is the bosonized 
form of the complex Weyl fermion ^'c that we constructed out of the two Major ana- Weyl edge modes in section UlI Al 
As we noted there, naturally couples to the (72 component of A. 

Gauge invariance tells us that the edge-modes must exist, but it does not determine their dynamics. We can, 
however, add a manifestly gauge invariant boundary term that ensures that the x field propagates unidirectionalUy 
at speed c = —A/kj. This term is 

J dxdt tr {[W~\d^+ A^)W][W-\d^ + c-^dt+ A^ + c-^At)W]} . (80) 

Including it, and writing W = e*'^^^, the effective action becomes 

^^p"(d,Ax) = ^ciA)-^wiv) + ^ J dxdt {d^x[{c-'dt + d.M} 
1 

+ 2^ 



■^eff 



^ J dxdt {{cAl + Al)d,;x} - ^ J dxdttr {{cAx + At)Ax} ■ 
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The terms containing x now constitute the action for a chiral boson interacting with the appropriate chiral component 
of the gauge field. The entire expression is invariant under gauge transformations U which reduce to the form e*'^^'^ 
on the boundary, and so maintain d = y there. We note that the factor of 1/2 in the "level" k is compensated for by 
a factor 2 coming from tr (cr|) = 2, so as to give the correct scale for a chiral boson representing a Weyl fermion. 

Since the action of the U{\) and SU{2) gauge groups commute with one another, the complete gauge invariant 
effective action containing all the low-energy degrees of freedom is the sum 

5^1* (A I*, d, Ax) = 25J>r(A'i>) + ^eT(d,Ax)- (82) 
Here the "2" in front of S'"g™ comes from the two spin components. 

V. CONCLUSIONS 

We have investigated the gauge invariance of the low-energy effective action for a 2 -|- 1 dimensional chiral superfluid 
coupled to external gauge fields. When the order parameter completely breaks the gauge group, as in the case of 
Abelian particle-number symmetry, the effective action becomes manifestly gauge invariant as soon as we include 
the bulk Goldstone modes among the dynamical fields. When the gauge symmetry is not completely broken, as in 
the case of spin-rotation symmetry, we found that additional, non-Goldstone, degrees of freedom were required for 
manifest gauge invariance. These were identified as being the spin-up and spin-down Majorona-Weyl edge fermions, 
which could be combined to produce a current that soaks up the remaining gauge dependence. 

The consequences of the effective actions for spin and particle-number currents also differ. The former has a true 
spin-Hall effect, with a dissipationless spin current proportional the spin "electric" field. The latter, we have argued, 
has only a "mock" Hall effect, the induced current being proportional to the change in density, and not to the external 
field causing the change. 
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APPENDIX A: DIRAC EQUATION 

The twisted-mass Dirac equation that results from our two-dimensional problem is a standard illustration of the 
theory of fractional charge [23. We review it here so as to make clear the contribution of the extended scattering 
states to the boundary current. 

We consider the one-dimensional Dirac Hamiltonian 

where A is a constant. 

We will compute the extra particle number that is accumulated in the vicinity of a; = when cf) is discontinuous, 
jumping abruptly from (j) — (j)^ when x < Q to cj) — (pR when a; > 0. Suppose the eigenstates of H are Xn with energy 
En- The ground-state number density is 

{^^^{x))=Y,\Xn[x)\\ (A2) 

where the sum is over occupied states, i.e. those with En < 0. Because the sum of |xnP over all states is independent 
of (f) by completeness, we can equally well write 



(Al) 



(V'V(a;)) const. - ^ |Xn(a:)P, 

B„>0 



(A3) 
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and this form is slightly more convenient. We will show that 

Q = j {ip''i){x))dx = ^{(j>R-M, 0<{<pR-(t>L)<'!T, 

= -^{<t>R-<t'L)-l, W < {(l)R- (Pl) <2w. 

When (p is constant we have a continuum of positive and negative energy eigenfunctions 



.1,4' _ pi(T34>/'2 



E + k + A 
E-k + A 



Akx 



2^/E{E + A) 

where E{k) = ±Vk'^ + A^. With the discontinuity present, we will have scattering solutions 



(A4) 



(A5) 



"_R 'r-k,E 



a 



(oul) I (pa 



R 



X > 0. 



The function ^ must be continuous at a; = 0, and from this condition we obtain the S'-matrix relation 



(out) 
(out) 



t r 
r t 



(in) 
(in) 



(A6) 



where 



t{k,E) 
r{k,E) 



1 



cos($/2) - [lE/k) sin($/2)' 

i(A/fc)sin($/2) 
cos($/2) - {iE/k)sm{^/2y 



(A7) 



Here $ is shorthand for (f)^ — (pR- 

In addition to the continuum states there is also a single bound state 



OC e" 



Eo+iK + A 
Eo-tK + A 

Ef) -iK + A 
Eo+iK + A 



e"'^^, X > 0, 
s«^. X < 0. 



The bound-state energy Eq is also determined by continuity at x = 0, which requires 

£^0 = Acos($/2), 
K = Asin($/2). 

These formula are valid for < $ < 2?:, where k is positive, and extend outside that interval with period 2tt. At 
$ = there is no bound state. As $ increases, a bound state peels off the upper continuum. It passes through E = 
at <1> = TT, and merges with the lower continuum as $ reaches 2tt. If $ increases beyond 27r, the process repeats with 
another state peeling off the upper continuum. Thus each 27r twist in $ results in the net transfer of one state from 
the upper continuum to the lower. 

Using the relation between Eq and k, we can write the normalized bound state as 



gi(0L+0H)/4 

-i(0i+0K)/4 



A complete set of states comprises xo together with 

^i^tE+r{k,E)i,tl^, 
t{k,E)4-^, 



Xk,E 



Xk,E = 



t{k,E)i>t:E, 



X < 0, 
X > 0; 

X < 0, 



A; > 0. 



fc < 



(A8) 



(A9) 



(AlO) 



k,E' 
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zA 

/Asin<I)/2 



■Rek 



-/AsinO/l 



FIG. 3: The contour F, showing the cut starting at ik = A. When < $ < tt there is pole at k — iAsin<l>/2, and a pole on 
the second sheet at fc = — iAsin<l>/2. 



These basis states therefore switch from a wave incident from the left to one incident from the right as k changes sign. 
After using |rp + |ip = 1, and the explicit form of the free eigenfunctions, we find 

E \xe{^)\'- f^^\XEA^)\'= const. + |^gr(fc,i?)(^^)e^»^N. (All) 

Here the "const." refers to terms that are independent of (f>R,L- We can improve the numerical convergence of the 
integral by using Jordan's lemma to push the contour of integration F into the upper half plane, as shown in figure 

We observe that the reflection coefficient has a cut running from A; = iA to ioo, and that for < $ < tt there is a 
pole in the upper half plane at A: = iAsin$/2. At ^ — tt the pole merges with the cut. For tt < ^ < 2tt the pole is 
apparently returning towards the real axis again, but a more careful investigation shows that it is now on the second 
sheet, and no longer contributes. (At the same time a pole at — iAsin$/2 has emerged onto the first sheet in the 
lower half plane. This is below the real axis, however, and also does not contribute.) 

The integral in ljAll|) then becomes 

- Asin<i>/2e-2'^l^-*/2|l-l + T , ^'/^"f, , ^ , ^ ^e^^-'-l. (A12) 

' 27T k2_A2 sin2($/2) Vk2-A2 ^ ' 

The first term, the pole contribution, is only present if < $ < tt. 

There does not seem to be a closed- form expression for the integral in ljA12p . but if we first integrate over x to get 
the total charge, we end up with an elementary integral 

dk A^sin^ 1 $ ^ 

^ = , -TT < $ < TT. (A13) 

A 27r k2 _ sin2($/2) - A2 2tt' ^ ' 

This expression extends to a 2tt periodic function of <f>, and the integral is therefore discontinuous at odd multiples 
of TT. After wc include the pole contribution, which is discontinuous at all multiples of tt, wc find that the total 
continuum contribution is discontinuous only at $ — (mod 2tt). Thus 

/°° dk <I> 
— \XE.k{x)\^ ^ const. + 0<$<27r. (A14) 

-oo 

The discontinuity at $ = (mod 2tt) is due to the sudden loss of the bound state from the upper continuum. As the 
bound state makes its way to the lower continuum, the spectral weight in the upper continuum is gradually recovered. 
We can now apply these results to compute 

j \Mx)?dx (A15) 



We observe that our sum over the positive continuum is equal to minus the sum over the negative continuum together 
with the bound state. The bound state has negative energy, however, if tt < $ < 27r (mod 27r). If we are computing 
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the ground state charge, and if < $ < tt (mod 27r), we should reduce the sum by unity. We also note that the 
constant in (|A14|I is fixed by the requirement that the the accumulated charge be zero when cj)]^ = (j)^. Thus 

< (0_R - (j)L) < TT, 

4>l) - 1, TT < ((Aii - (/-l) < 2tt. (A16) 

This result repeats with 2tt periodicity. 

This expression for Q is consistent with the basic result of Goldstone and Wilczek , that a total charge of 




(A17) 



is drawn into a region as the phase (f> is slowly twisted. The difference between (|A16ll and (|A17ll arises because the 
latter does not keep track of individual particles that are lost to the reservoir when their energy exceeds the chemical 
potential. 

In a superfluid the role of charge and current is interchanged. The fractional charge, multiplied by fc/ and divided 
by two, gives the momentum density, or equivalently, the mass current psj . 



Q = 



1 

1 

2^' 
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